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Abstract 

Wc define a new type of self-similarity for one-parameter families of stochastic processes, 
which applies to a number of important families of processes that are not self-similar in the 
conventional sense. This includes a new class of fractional Hougaard motions defined as 
moving averages of Hougaard Levy process, as well as some well-known families of Hougaard 
Levy processes such as the Poisson processes, Brownian motions with drift, and the inverse 
Gaussian processes. Such families have many properties in common with ordinary self-similar 
processes, including the form of their covariance functions, and the fact that they appear as 
limits in a Lamperti-type limit theorem for families of stochastic processes. 

Key words: Exponential tilting; fractional Hougaard motion; Hougaard Levy process; Lam- 
perti transformation; power variance function. 

Mathematics Subject Classification (2010): Primary 60G18, 60G22; Secondary 60F05 

1 Introduction 

The purpose of this paper is to extend the notion of self-similarity to one-parameter families 
of stochastic processes. Recall that a real-valued stochastic process {X(t),t > 0} is called self- 
similar with Hurst exponent H £ (0, 1) if it satisfies the scaling property 



X(ct) = c M X{t) for t > 0, (1.1) 



for al l c > 0, where = denotes equality of the finite-dimensional distributions ( Embrechts and Maejima 



2002 ) . The corresponding process with drift /i€K, defined by 

X(p;t)=X(t) + /jtfart>0, (1.2) 

does not, however, satisfy (jl.ip for \x ^ 0, an example being Brownian motion with drift. Instead, 
the family of processes (|1.2p satisfies the following scaling property, 

X(nc H - l ;cb) = c H X{fi; t) for t > 0, (1.3) 
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for all c > 0, and we shall hence propose (jl.3p as a new definition of self-similarity for families 
of stochastic processes, cf. Definition 12.11 below. 

We show that there are in fact many important families of stochastic processes that satisfy 
(|1,3|) without having the drift form (\1.2\) , and for values of H that d o not necessarily belong t o 
(0, 1). One such example is the class of Hougaard Levy processes ( Lee and Whitmore . 19931 ) . 
which includes for example the family of Poisson processes {H = 0), certain gamma co mpound 
Poisson processes {H < 0), and the family of inverse Gaussian processes (H = 2), cf. IWasanl 
(1968). A further example is a new class of fractional Hougaard motions defined as moving 
averages of Hougaard Levy processes, generalizing fractional Brownian motion. As we shall see, 
such processes have many properties in common with ordinary self-similar processes, as reflected 
for example in the familiar form of their covariance functions. This represents an important step 
forward compared with conventional self-similar processes, where the only processes with finite 
variance are the fractional Brownian motions. 

In Section [21 we present the new notion of self-similarity for families of stochastic processes. 
We show that their covariance structure is completely determined by the so-called variance func- 
tion, and we study the role of power variance functions. In Section [3] we explore the connection 
between self-similarity and exponential tilting, and show the self-similarity of Hougaard Levy 
processes. In Section H] we show the self-similarity of the class of fractional Hougaard motions 
using their moving average representation. In Section [5l we introduce a Lamperti-type transfor- 
mation, which transforms a self-similar family into a family of stationary processes. In Section [6] 
we consider some Lamperti-type limit theorems, where families of self-similar processes appear 
as limits of suitably scaled families of stochastic processes. Finally, in Section [7J we investigate 
the case H = 1 and its relation with exponential variance functions. 



2 Self-similarity 
2.1 Definition 

We consider a family of real- valued stochastic processes X = {X(fj,; t) : fj, G O, t > 0} , all defined 
on the same probability space, where t denotes time, and where the family is indexed by a real 
parameter \x G £1, where Q is an interval satisfying R + C £1 and cfl C for all c > 0. We allow 
0, to contain the values ±oo, and define = intS7, the interior of f2. These conventions are 
important for the examples that we will discuss. Motivated by the above discussion of processes 
with drift, we now propose an extended definition of self-similarity for families of stochastic 
processes. In the following we use the shorthand notation X{n;-) = {X(fi;t) : t > 0} etc. to 
denote a particular stochastic process in X. 

Definition 2.1 (i) A family of stochastic processes X is called self-similar with Hurst exponent 
H G R and rate parameter fi G f2 (H-SS) if 

X^c 11 ' 1 ; ct) = c H X(fi; t) for t > 0, (2.1) 

for all c > and fi e !1. (it) We say that X has stationary increments if each process in the 
family has stationary increments, that is, for all s > and [x G £l 

X{fi; s + t)- X(n; s) = X{fx; t) - X(fi; 0) for t > 0. (2.2) 
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The family X is called H-SSSI if it is H-SS and has stationary increments. 



When necessary, we refer to the definition (|2.ip as general self-similarity, whereas the con- 
ventional definition ([Lip for stochastic processes will be called strict self- similarity from now 
on; a terminology borrowed from stable distributions. The motivation for expanding the use of 
the term self-similar in this way is that many properties of strictly self-similar processes have 
immediate analogues for general self-similar families, in particular for the covariance structure 
(cf. Section EZD- 

A family of processes with drift X = {X(t) + fit : [i £ R, t > 0} is general sense H-SS if and 
only if the stochastic process X(-) is strictly self-similar with Hurst exponent H. It is clear that 
general self-similarity constitutes an extension of strict self-similarity, since if X satisfies ([2. 1 1) . 
then the distribution of X(fx; t) — fit will generally depend on fj,, as illustrated by the examples 
in Sections 13.21 and HI Nevertheless fi plays the role of rate for the process X(p,;-), in the sense 
that the dimension of \x is the unit of the process X(fi; •) per unit of time, thereby providing an 
extension of the idea of a drift parameter. When applicable, the values \x = and \x = ±oo in 
(|2.ip correspond to the strictly self-similar processes X(0; •) and X(±oo; •), respectively. 

An equivalent way of writing the definition (|2.ip is as follows, 



C - H X(fic H - 1 ;ct) = X(ji; t) for t > 0, (2.3) 

for all c > and \i € O. Noting that the right-hand side of (|2.3p is independent of c, we observe 
a collapse of the scaled marginal distributions of X on to a single distribution independe n t of c , 
much like the notion of universality in turbulence, see Barndorff-Nielsen and Schmiegel ( 20081 ) 



and references therein. Taking c = 1/t in (|2.3p . we obtain the useful relation 

^X^t 1 ' 11 -!) ~ X((i;t) for t > 0, (2.4) 

for each fj, E Q, where ~ denotes equality of the marginal distributions for each t. Given Hgt, 
the left-hand side of (|2,3p represents a multiplicative transformation group {R? : c> 0} , defined 
by 

R?X(p; t) = c- H X{juF-\ct). (2.5) 
This is a continuous analogue of t he ren ormalization group in the sense of Jona-Lasiniol ( 19751 ). 



see also Embrechts and Maeiima ( 20021 . p. 15), with fixed point given by (j2.3|) . a topic that we 



shall return to in Section [6l 

The general interpretation of the definition ([ 2 . 1 1) is that a rescaling of X(fi;t) is equivalent 
in distribution to a simultaneous rescaling of time t and rate fj,, so that non-trivial solutions 
to (|2.ip require that the different processes of the family are suitably linked, for example by a 
drift term, as in (|1.2p . or by means of exponential tilting, cf. Section [3l In fact, (|2.4p suggests 
that a simple rescaling t H X^t 1 ^ 11 ) of a given process X(-) provides a trivial solution to ()2.ip . 
Similarly, for H ^ 1 and /i/0we may take c = \fi\ l ^ 1 ~ H ^ 1 in (|2,3p . which yields 

X(fi; t) = \fx\ H/{H - t] X(sgn(/i); H V(1 ~ H) *) for t > 0, 

where sgn(^i) is the sign of \x. Hence, if X(-) is a given stochastic process, then the rescaled 
family defined by X(p,;t) = fj, H ^ H -^X(/j}^ 1 ~ H h) for /z,t > is also a trivial solution to (|2TTj) . 
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Note that taking t = in (pTTJ) gives X{fxc H ~ 1 ; 0) ~ c H X(fi; 0), which shows that if 0) ~ 
for one value of fx, then the same is the case for all values of fx with the same sign. In most 
cases we have X{fi; 0) = 0. 

2.2 Covariance structure 

Let X be an .ff-SSSI family of stochastic processes. If X(fx; 1) has finite second moments for all 
fi £ Cl (second-moment assumptions) then fj2.4[) implies that X(fx;t) has finite second moments 
for all t > and fx E Cl. We shall now explore the covariance structure for such families. Most 
of the results in the following require H ^ 1, and the case H = 1 will be considered separately 
in Section [71 The next result concerns the structure of the first moment. 

Proposition 2.1 Let X be an H-SSSI family of stochastic processes with H ^ 1, and assume 
that X(fx;t) has finite expectation for all fx 6 Cl and t > 0. // the function t t— > E[X(//;t)] is 
continuous on [0, oo) for each fx E Cl, then there exist functions a and b such that 

E [X(fx; t)] = a(sgn(») \fi\ H/{H - l) + 6(sgn(/x)) \(J,\ t for t > and fx G U. (2.6) 

If CI = M and < i? < 1, £/&en continuity of the function \i i— >• E [X(/i; £)] implies a = 0. 

Proof: For given fj, E. CI, the assumption of stationary increments (|2.2[) implies that, for s, i > 0, 
E [X( M ; s + i)] + E 0)] = E [X(p; s)) + E [X(/x; t)] . 



Using standard results for the Cauchy functional equation (cf. iBingham et al.l . 119871 . p. 4), the 
continuity of E [X(fx; •)] implies that 

E [X(p; t)] = a(/i) + b(p)t for t > 0, (2.7) 

for suitable function a and b defined on Cl. In view of (|2.1|) . we obtain, for all c, t > and fx E. CI, 

a{Lic H ' 1 ) + b{Lic H ' l )ct = c H a{Li) + c H b(Li)t, 

which implies that 

a(fic H - 1 ) = c H a(fi) and bine 11 " 1 ) = c H_1 6(/i). 

By inserting ii = ±1 and x = ±c^ _1 in these equations, we obtain a(x) = \x\ H ^ H ~^ a(sgn(x)) 
and b{x) = \x\ 6(sgn(x)), in agreement with (|2.6|) for fx ^ 0. 

If Cl = K, such that G Cl, then we have already seen in Section f2.il that the process X(0; ■) 
is strictly self-similar, in which case it is well-known that ^(0; 0) = 0. This is in agreement with 
the limit of (|2.6p as fx — > 0, provided that either a = or \fx\ H ^ H ~^ — > as fi — > 0, the latter 
being the case outside the interval < H < 1. Hence in the case < H < 1, continuity of the 
function /i i— >■ E [X(//;t)] implies a = 0. This completes the proof. ■ 

The constant term a(sgn(//)) may be removed from the processes by subtraction. 

Thus, if the family X satisfies the conditions of Proposition 12.11 we define the corresponding 
centered family Xq by 

X (fi; t) = X(fi; t) - a(sgnOu)) \fx\ H ^ H ^ . (2.8) 
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It is easy to check that the family Xq is again ii-SSSI, now with a(sgn(/i)) = 0. In many cases 
the mean function (|2.6p for the centered family Xq takes the following simple form 

E [X (fi; t)] = /it for t > and /i£fl, (2.9) 

corresponding to b(fi) = fi. The latter form may be assumed without loss of generality (i.e. up 
to a rescaling of ji) in the important case Vt = R + whereas for J) = R this requires 6(1) = b(— 1). 
From now on we assume that the family X is centered, unless otherwise stated. 

We shall now derive the covariance structure of the family X under second-m oment assump - 
tions, which is done in much the same way as for strictly self-similar processes (cf. iTaqqul . 12003) . 
Let us introduce the variance function V : fi — > R + defined by 

V(fJ,) = Var [X(n; 1)] for fi € 0. (2.10) 

Consider a given value of fi G f2. The marginal variance of the process may be obtained by 
calculating the variance on both sides of (|2.4p . giving 

Var [X(ji; t)} = t 2H V{fit l ~ H ) = V H (fx; t), (2.11) 

say, for all i > 0. The whole covariance structure of X may now be expressed in terms of the 
variance function V. Thus, for s, t > we obtain 

Var [X(fi; t) - X(n; s)] = Var [X(fi; s)] + Var [X(p; t)] - 2Cov [X(fi] s),X(fi; t)] . (2.12) 

Using the stationarity of the increments together with f)2.11f) . this gives the following covariance 
function for the process X{pL\ ■) for s,t > 0, 

Cov [X{n; s), X(ji; t)) = \ [V H (//; s) + V H (/x; t) - V H ( M ; \t - s\)] . (2.13) 
2.3 Power variance functions 

Let us consider the important case where a centered -fT-SSSI family of stochastic processes X 
has power variance function V((J.) = a 2 /jP for some p € R and a 2 > 0, where p is called the 
power parameter. Here we take the interior of the domain for [i to be Oo = R for p = and 
Q p = R+ for p ^ 0. In this case, the variance (|2,lip becomes 

Var [X(ji; t)] = a 2 fi p t 2H+p ^- H ^ = a 2 fi p t 2 ~ D for t > 0, (2.14) 

say, where D is the fractal dimension defined by 

D = (H-l)(p-2), (2.15) 

with domain D £ [0,2] (cf. Lemma 12.21 below). 

The covariance function (|2.13p now takes the form 

Cov[X(fJ,;s),X(n;t)] = a 2 fi p R D (s,t) for s,t > 0, (2.16) 
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where 



R D (s,t) 



1 



S 2-D + t 2-D 



\2-D 



for s, t > 0. 



Compared with (|2,13p . the covariance function (|2.16|) is now a product of the variance function 
a 2 \JP and a function depending on s and t only. Let us also consider the correlation between two 
non-overlapping increments, which may be expressed in terms of r = \fs]i, 



Corr [X{fi] s),X(ji) s + t)- s)] 



1 r 



- X + r) 



2-£> 



„2-D 



-2+D 



fors,t>0. (2.17) 



Lemma 2.2 T/ie function Rd is non-negative definite if and only if < D < 2. 



Proof: Let us first take p = in (|2.15p . corresponding to D = 2 (1 — H) , in which case 



#2(i-ff)(M) = g 



,2H 



27^ 



i 



,2H 



has the same functional form as the covariance function o f standard fractional Brownian motion , 
which is known to be non- negative definite for < H < 1 ( Tag qui . 20031 : lEmbrechts and MaeT 



uma 



2002). It follows that Rd is non-negative definite for any D € [0,2]. Let us now show the ne- 
cessity of this condition. Since the correlation (|2.17p must be less than or equal to 1, we obtain 
for r = 1 that (2 2 ~ D — 1 — l) /2 < 1, which is equivalent to D > 0. Now consider the case 



D > 2, where the function on the right-hand side of (|2.17p behaves asymptotically as 



-2+D 



/2 



as r — > oo. For r large this is in contradiction with (|2.17p being a correlation, so that we must 
have D < 2. ■ 



Lemma 2.3 Let X be an H-SSSI family with power variance function. Then 

1. The correlation {2.11 ) is positive for < D < 1 and negative for 1 < D < 2. 

2. The process X([a; •) has uncorrelated increments if and only if D = 1, where A2.16\) becomes 

Cov [X(ji; s),X{n; t)} = ^o-V mm {s, t} . (2.18) 



Proof: Item[TJ The correlation (|2.17p is positive if and only if 

(r- l +rf- D >r 2 - D +r- 2+D . 

In terms of the strictly convex function f(x) = log {e x + e~ x ) , this is equivalent to the inequality 

(2-D)f(x) >f(x(2-D)). (2.19) 

Using the convexity of / along with the fact that f{x) behaves asymptotically as \x\ for \x\ large, 
we find that the inequality (|2.19p is satisfied for < D < 1, whereas the opposite inequality 
holds for 1 < D < 2, which proves Item[TJ 

Item[2j The correlation (|2.17p is zero if D = 1, which by Item [1] is also necessary. ■ 
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V 


Positive correlation 


Negative correlation 


V < 2 (a i [0, 1]) 
p>2 (a €(0,1)) 


l/a<H<\ 
1 < H < 1/a 


(2 - a) /a < H < 1/a 
1/a < H < (2- a) /a 



Table 1: the sign of the correlation (|2.17p as a function of a and H. 



It is convenient at this point to introduce the parameter a € [—00,00), defined by the 
following one-to-one transformation of p £ (—00, 00], 

a = a(p) = l + (l-p)-\ (2.20) 

with the conventions that a(l) = —00 and a(oo) = 1. By Item [2] of Lemma 12.31 we find that 
the correlation between two non-overlapping increments (|2.17p is zero if and only if H = 1/a, 
with the convention that a = — 00 corresponds to H = 0. A further reason for our interest in 
the parameter a comes from the connection with a-stable distributions in the case a G (0,2], 
see Section 13.21 where we discuss a class of self-similar families of Levy processes. 

We may now express the domain for the Hurst exponent H, for given q, as an interval with 
endpoints 1 and (2 — a) /a, which follows from the domain < D < 2 via (|2.15p and (|2.20p . 
This domain is summarized in Table [H along with the sign of the correlation (|2.17p (Item [T] of 
Lemma f2.2p . By comparison, in the case of strict self-similarity, the Hurst parameter is restricted 
to the domain H > 0, and to the smaller domain < H < 1 under second-moment assumptions, 
cf. iTaaaJ tooi ). 



The importance of power variance functions will become clear in Sections 13.21 and HI where 
we give examples of processes with independent increments and with correlated increments, 
respectively, giving rise to .ff-SSSI families corresponding to many of the possible combinations 
of H and p discussed above. We note in this connection that D = and H ^ 1 implies p = 2 
(a = 0), in which case the covariance function (|2.16p takes the form 



Cov [X (fj,; s),X(p; t)] = a 2 fi 2 st for s, t > 0. (2.21) 



We shall return to this case in Section [3.2i The other case where D = 0, namely H = 1, will be 
considered in Section [71 



3 Exponential tilting and Hougaard Levy processes 



3.1 Exponential tilting 

We now consider t he important case wher e X is a natural exponential family of centered stochas- 

Let P* denote the marginal distribution of 
and denote expectation under P M 
and P* , respectively. Let k denote the cumulant generating function of , defined by 



tic processes, cf. iKiichler and Sorensenl (|1997I ). 
X(/i; t) under the probability measure P^, say, and let E 



k(0) = log Ei 



(3.1) 
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with effective domain = {6 £ R : k(9) < oo}, assumed to have non-empty interior. Assuming 
that P t = {P* : \i £ f2} is a naturai exponentiai family for t = 1, we shall see that this is then 
the case for all t > 0. We then call X an exponential PT-SSSI family. 

Let us consider the exponential change of measure from Pi to P«, corresponding to the 
likelihood ratio 

K„{t) = exp [9\X(l;t) - k\(0\)] , (3.2) 

where is the cumulant generating function for P* and the canonical parameter 9\ is chosen 
such that Ei [X(l;i)] = t. The exponentially tilted probability measure P^ may then be defined 
from Pi as follows: 

Pf l (A)=E 1 [l A K ll (t)], (3.3) 
where 1a denotes the indicator function for the Borel set A. More generally we have 



E< [tf]=Ei [UK^t)], 



(3.4) 



where U is a non-negative random variable. 

We shall now derive the cumulant generating function for general \i and t = 1. Standard 
exponential family arguments along with (|3.4p yield 



kJ(u) 



«X( W 1) 



with effective domain 
domain f2 = k(0), and 



logE M 
= log Ei 

= log Ei r e (^+«)x(i ; i)- K (el) 

= K (u + e\)-K(e\), (3.5) 

?}. The mean and variance function of X(/j,;l) are \i = k{9\) with 



V(p) = ko k' 1 ^) for (j, en, (3.6) 

respectively, where dots denote derivatives. The assumption that P/ has mean 1 implies that 
k(0) = 1. It is well known that t he variance funct ion V together with its domain $7 characterize 
the exponential family (j3.5|) . cf. IJorgensenl (jl997l . Ch. 2). 

Turning now to the case of general t > 0, we obtain by means of (|2.4p the following expression 
for the cumulant generating function «* 



logE^ 

logE„ [exp (u^A-(Mt 1_H ;l))] 

/t H ))-K{t H {el/t H )) 



k (t H (u + 0*, 



(3.7) 

/i. In particular, 

this implies that X(n;t) has mean fit and variance given by f)2.11f) . in agreement with the results 
of Section I 



for values of u such that ut H + 9 f £ 0, where 9 l is the solution to k(9 u )t 



H-i 
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The form (|3.7p shows that P is indeed a natural exponential family for each t > 0, corre- 
sponding to the cumulant generator n{t H •) and canonical parameter 9^/t H . We may interpret 
(|3.7p as representing a translation structure in the Fourier domain, as compared with the drift 
term /it in (|1.2|) . which represents a translation structure in the sample space. 



3.2 Self-similarity of Hougaard Levy processes 



We now consider the case where X is a centered ff-SSSI family of Levy processes. Under 
second-moment assumptions, the independence of the increments implies the following form for 
the covariance function: 



Cov [X(/j,;t),X (fx; t + s)] = Vax[X(ji;t)] = tV(jj) for s,t > 0. 



(3. 



Under the further assumption that the family X is exponential, it is clear from the results of 
Section that the variance function V with its domain f2 characterize the family X among 
all exponential f/-SSSI families, because V characterizes the exponential family of marginal 
distributions of X(n; 1), which, in turn, characterizes the corresponding family of Levy processes. 
From now on, we denote the marginal distribution of t~ l X (fi\t) by the sym bol ED(/x,t), 



where ED(//, t) is an exponential dispersion model in the sense of |j0rgensenl (|1997l . Ch. 3). This 
distribution has mean fi and variance t -1 U(/x) for t > and /j£fl. The above characterization 
of the family X in terms of the variance function V may now be rephrased as follows. Namely, 
the variance function V with domain 0, char acterizes th e expo nential dispersion model ED(/i, t) 
among all ex ponential d i spersi on m odels, cf. J0rgensenl (1997, Ch. 3). 



Following I J0rgensenl (|1992l ) and lLee and Whitmord (|1993h . we define a Hougaard Levy pro 



cess S p (jj,]t) to be an exponent i al fam il y of Levy process with power variance function V (n) = 
a 2 n p , see also iHougaard et al.1 (|l997l ). |j0rgensenl (|l997l . Ch. 4) and IVinogradovl (j2Q0j ). We 
define the Tweedie distribution Tw p (/i,t) to be the corresponding exponential dispersion model, 
such that 

t^Spifi-t) ~Tw p (/x,t), (3.9) 

with mean fi and variance t~ x o 2 pP (note that we suppress the parameter a 1 in the notation). 
The parameter domains are p £ A 



-- M\(0, 1) 


a 2 > 0, and \x G O p , 


' [0,oo) 


for p < 


R 


for p = 


R + 


for 1 < p < 2 


, (0,oo] 


for p > 2. 



The Tweedie distribution satisfies the following scaling property for all p G A, /i G fi p and 
i > 0, 

cTw p (/i, t) = Tw p (c/i, c p ~ 2 t) for c> 0, (3.10) 

which characterizes the Tweedie model among all exponential dispersion models, see Jorgensenl 
(j 19971 . Ch. 4). The Hougaard Levy process S p (fi;t) is hence characterized by the following 
distribution of the increments, 



S p (/j,;t) ~ tTw p (fi,t) = Tw p (fit,t p 



(3.11) 
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with variance ta 2 fi p . We shall now characterize the Hougaard Levy processes in terms of self- 
similarity. Recall that the parameter a is defined as a function of p by (|2.2(jp . and note that 
p € A corresponds to the domain a G [—00, 1) U (1, 2]. 

Theorem 3.1 The family of Hougaard Levy processes S p (p>; •) may be characterized as follows. 
Let X be an exponential family of Levy processes. Then X is self-similar with Hurst exponent 

H=-, (3.12) 

a 

with the convention that a = —00 (p = 1) corresponds to H = 0, if and only if X has power 
variance function with power p 7^ 2. 

Note that the value H = 1/a is in agreement with the condition for uncorrelated increments 
of Lemma E21 Item[2j 

Proof: Let us first show that the Hougaard Levy process S p (/j,;t) is 1/a-SS. Since s Hougaard 
Levy process has stationary and independent increments, it suffices to consider the marginal 
distribution of S p (fi; t), so that we must show 

S p {nc H ~ l ]ct) ~ c H S p {fi;t) for t > 
with H = 1/a. Using (|3. 1 1 1) . the left-hand side has marginal distribution 

S p {nc H - l ;ct) ~ Tw P 0utc H , (ct)^ 1 ). 
The scaling property (|3.10p yields the following marginal distribution for the right-hand side, 

c H S p (fi;t) ~ c II Tw p (pkt,t p - 1 ) 

= 1vj p {^itc H ,c H ^- 2 hP- 1 ) 

= Twp^,^)"- 1 ), 

where we have used the fact that H (p — 2) = p — 1, which follows from (|2.20p and (|3.12p . This 
implies that S p (n; t) is iT-SSSI with H given by f|3. 12j) . 

To show the reverse implication, let us assume that the exponential family of Levy pro- 
cesses X is self-similar with Hurst exponent H = 1/a. By (|3.8p the variance of the process is 
Var [X(n; t)} = tV(fj,), and comparing with (|2.1ip . we obtain the equation 

tV(n) = t 2H V{nt l ~ H ). (3.13) 

By taking \i = 1 in f|3. 13|) and redefining [i to be [i = t 1_H , we obtain the following solution: 

v(v) = v(i)^ 1 - 2H ^ 1 - H \ 

and using (|3.12|) and (|2.20p we find the power to be (1 — 2H) / (1 — H) = p. Hence V is a power 
variance function with power parameter p and a 2 = V(l). ■ 
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4 Fractional Hougaard motion 



We shall now consider stochastic integration with respect to a Levy process, which in turn will 
be used to construct a class of fractional Hougaard motions as moving av erages of Hougaar d 
Levy processes. Some of the results in the following are parallel to results of Marquardt ( 20061 ) . 
but unlike Marquardt, who considers Levy processes with zero mean, finite variance, and no 
Brownian component, we shall consider stochastic integration with respect to an arbitrary Levy 
process. 



4.1 Stochastic integration with respect to a Levy process 

W e now consider stochastic integration with re spect to a Levy process X. We follow the approach 
of Barndorff-Nielsen and Thorbi0rnsen ( 20061 ) , which has the advantage that we work directly 
with the cumulant transform (log characteristic function) and the cumulants, although there are 
mor e general approaches t o stoc hastic integration with respect to a Levy process available, see 
e.g. Rajput and Rosinski ( 19891 ). 

We shall define the stochastic integral 



f(u)dX(u), 



(4.1) 



for a non-random real function /. Initially we consider the case where < a < b < oo, 
after which we extend the integration interval to the whole real line. We denote the cumulant 
transform of an infinitely divisible random variable Y by 

C{z;Y) = log E[e izY ] . 

Due to the infinite divisibility, the characteristic function E [e 12 ^] has no zeroes, and hence C 
is well-defined by means of the principal branch of the complex logarithm. Recall that 

E [Y] = -iC{0; Y) and Var [Y] = -(7(0; Y), (4.2) 

provided that C is twice differentiable at zero, where C(z;Y) and C(z;Y) denote the first and 
second derivatives of C, respectively, with respect to z. 

Lemma 4.1 Assume that < a < b < oo. Let X be a Levy process, and let the function 
f : [a, b] — > R be continuous. Then the stochastic integral ft4-l\ ) exists in the limit, in probability, 
of approximating Riemann sums. The distribution of the random variable j^. 1\ ) is infinitely 
divisible, and has cumulant transform given by 



C[z]J f(u)dX(u)j=J C(zf(u);X(l)) du for z G R. (4.3) 
If X(t) has finite second moments, the mean and variance of the random variable j^. 1\ ) are 

E I Hu)dX(,i) = E[X(1)] [ f(u)du, (4.4) 



Var 



f{u) dX(u) 
f(u) dX{u) 



Var[X(l)] / f 2 {u)du. 



(4.5) 
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Proof: The stochastic integral (14. ID , and the result (|4.3p follow from Lemma 2.5 of 
Barndorff-Nielsen and ThorbjOrnsen ( 20061 ). By differentiating (|4.3p and putting z = 
obtain 



wc 



E 



f(u)dX{u) 



[ b f(u)C(0; 

J a 



X(l)) du 



E [*(!)] / /(«)d«, 



where we have used the fact that — iC (0;X(1)) = E[X(1)]. Differentiating (|4.3p twice and 
putting z = we obtain 



Var 



/(u)dX(u) 



/ 2 (n)(7(0;X(l)) d« 
VarLY(l)] / f 2 (u)du, 



where we have used the fact that —C(0;X(1)) = Var [X(l)]. This shows the last two results. In 
particular, the two integrals on the right-hand side of (|4.4p and (j4.5[) are finite, each being the 
integral of a continuous function on a compact interval. ■ 

We shall now extend the integration interval to the whole real line. To this end we first need 
to extend the time domain of the process X to R. Since we are dealing with a Levy process, we 
have that X(0) = 0. We then extend the process to the negative half-axis by assuming 



X(-t) = X(0) - X(-t) = X(t) for t > 0, 



(4.6) 



cf. iTaaaul l|200ah 



Proposition 4.2 Assume that —oo < a < b < oo. Let X(t) be a Levy process on M., and let 
f : (a, b) — > R be continuous on (a,b) \ {0}. Assume that 

VzeR: / \C(zf(u);X(l))\ du < oo. (4.7) 

J a 

Then the stochastic integral J^.i[ ) exists in the limit, in probability, of the sequence 

f(u)dX(u)\ , 

where a n and b n are arbitrary sequences in (a, b) such that a n < b n for all n and a n \ a and 
b n f b as n — >• oo. The distribution of the random variable |^.i[ ) is infinitely divisible, and 
has cumulant transform given by |J3p. If X(t) has finite second moments, then the mean and 
variance of the random variable (CO) are given by \4-4\) an d H-^\)> respectively. 
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Proof: First assume that < a < b < oo. In this case, the stochastic in t egral and the result 
follow from Proposition 2.6 of Barndorff- Nielsen and Thorbiornsen ( 20061 ). In the case 



-oo < a < b < we make the substitution v = —u, and rewrite the integral as follows: 

f(u)dX(u)= [ b f(-v)dX(-v)= [ a f(-v)d[-X(-v)}. (4i 



b 



In view of (|4.6p . the process {— X(— t) : t > 0} is a Levy process on the positive half-line. Since 
< — b < —a < oo, the right-most integral of (|4.8p is defined according to the first case just 
considered, thereby providing the required definition of the left-most integral of (|4.8p . In the 
case where a is negative and b is positive, we split the integral in two parts, 



/ 

■J a 



b f(u)dX(u)= f°f{ U )dX(u)+ f b f(u)dX(u), 

J a JO 



and use the fact that both integrals on the right-hand side of the equation are now properly 
defined. The remainder of the results follows from Lemma 14.11 bv elementary arguments. ■ 



4.2 Fractional Hougaard motion 



We shall now define a class of fractional Hougaard motions by means of stochastic integration 
with respect to a Hougaard Levy process. This approach is similar to the moving average 
representation of f ractional Brownian motion , wher e the integration is with respect to ordinary 
Brownian motion ( Mandelbrot and van Nessl 19681 ) . Fractional Hougaard motions provide our 
main exampl es of npn-Ley y self-similar fa milies o f stoch astic processes. 

Following lieTIr] (|l994L pp. 56-59) and lTaaaul (|2003h . we define the weight function Wh(t, u) 
for t > and uGtby 



w h (t,u) 



u) h - {-u) 
(t - u) h 




for 
for 
for 



u < 
< u < t 
t < u, 



(4.9) 



where h is a real constant. Using the notation u + 
function in the more compact form 



max {0, u} we may also express the weight 



w h (t,u) = (t-u)+ - (-«)- 



(4.10) 



For a given p £ A, we let {S p (/j,;t) : t G R} denote the Hougaard Levy process defined in 
Section 13.21 extended to the whole real line by means of 
and we recall from (|2.20p that a = 1 + (1 — p)^ 1 . 



For simplicity we take a = 1, 



Proposition 4.3 Assume that p G A\{1,2} and H < 1/a, and define the fractional Hougaard 
motion 5 Pi ij(/x;f) by the stochastic integral 



uih(t, u) dS p (fj,; u) for t > 0, 



(4.11) 
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where h = H — 1/a and fj, E £l p . Then the family S p ^h{h] t) is H-SSSI with cumulant transform 



n — 1 roc 



a 



1 + 



izw h (t,u) 
{a - l)//V(<*-i) 



1 \ du for z G M. (4.12) 



For 1/a — 1 < H < 1/a, the process Sp^i^t) has mean zero, and for 1/a — 1/2 < H < 1/a, 
the process has variance 



/OO 
-oo 



v) dv. 



(4.13) 



Proof: Our method of proof is inspired by iBeranl ([1994], pp. 56-59), and we shall hence show 
self-similarity by a direct argument, rather than via the covariance function, say, as is usually 
done for fractional Brownian motion. We first need to check that the integral fj4.7|) is finite for 
the function f(u) = Wh(t,u). Let us define the complex analytic function K a by 



K a (x) 



a — 1 



a 



a — 1 



for Rc ■ 



a — 1 



> 



(4.14) 



(jjprgensenl . 119971 . p. 131), where the domain may be extended to the imaginary axis when a > 0, 
and to the whole complex plane when a = 2. Then S P ((J,;1) has cumulant transform 



1 



iz 



1 



C (z; S p (n; 1)) = K a (6 + iz) - K a (9) = K a (9) 
where 6 = (a — l)i± l ^ a ~ l \ The integral ()4.7p hence takes the form 

/>oo 

\C (zwh{t, u); S p (n; 1)) | du = / \n a {9 + izwh{t,u)) — K a {0)\ du. 



In the limit u — > — oo, a Taylor expansion of K a yields 

\n a {6 + izw h (t,u)) - K a {6)\ = O (\w h (t,u)\) , 

which implies integr ability in this limit, due the assumption that h = H — 1/a < 0. For u 
near zero or i, the integrability of \K a (9 + izwh(t, u)) — K a (9)\ = O (\wh(t,u)\ a ) requires ah = 
aH — 1 > —1, or equivalently Ha > 0, which is satisfied because a < 2 and H < 1/a. By fj4.3j) 
we find that the process S Pt H(t) has cumulant transform 



C(z;S p , H {lJ>;t)) 



C (zw h (t,u);S p (fi;l)) du 
[n a (9 + izwh(t, u)) — K a (9)] du 



K a (9) 



IZ 



1 + -^w h {t,u) 



1 > du, 



which implies (|4,12p . in view of the above definition of 9. 
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Now let us calculate the mean and varia nce of SgjfU t'jt) using (|4.4p and (|4.5p . Using well- 
known results about the integral of Wh (e.g. Beran . 19941 . formula (2.22) p. 59) we obtain that 
the mean is finite for 1/a — 1 < H < 1/a and 

/DO 
w ft (t,iOdu = 0. 
-DO 

To calculate the variance we use the following scaling relation for the weight function ?%, 

w h {t,u) = t h w h (l,ut- 1 ). (4.15) 
Using the substitution v = ut~ l , and with h = H — 1/a, we obtain for 1/a — 1/2 < H < 1/a, 

POD 

Var [S p H (n;t)] = \i v / w\(t,u)du 



oo 

DO 



i/ — oo 

/■oo 



In order to show that the family of fractional Hougaard motions 5 p #(/i;i) is -ff-SS, we shall 
use the fact that the Hougaard Levy process S p (fj,; t) is 1/a-SS with a = 1 + (1 — p)^ 1 - Consider 
(|4,lip with arguments fxc^ -1 and ct. We note that the function Wh satisfies the following scaling 
relation similar to (|4.15p . 

w h (ct,u) = c h w h (t,uc _1 ) for c> 0. (4.16) 

In view of (|4.16|) . we obtain 

/oo 
w h {ct,u) d5 p (/ic H_1 ;u) 
-oo 

/oo 
w h (t, uc~ l ) dS p {nc H ~ 1 ;u). 
-oo 

Making the substitution v = uc~ l , this becomes 

/oo 
w h (t,v) dS p (ficP~ 1 ; cv) . 
-oo 

Recalling that h = H — 1/a, and using the 1/a-SS property of S p (/i;t), this implies 

, roo 

S p , H {p^-\ct) = c H - l l a / w h {t,v)c l / a dS p (fi;v) = c H S p>H (Kt). 

J — oo 

Hence, S p H(fj>;t) is self-similar with Hurst exponent H. 
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To show that SpHfat) has stationary increments, we observe, using the expression (I4.10p 
for Wh, that for any s, t > 0, 

/oo 
{w h (t + s,u)- w h (t, u)} dS p (ji; u) 
-oo 

{(t + s-u)%- (*-«)+} dS p (fi;u) 



oo 
oo 



oo 
oo 



r U* -V)% -(-V)%\ dS, fa V) 

S p ,H(fJ>;s), 



where, in the second-last equality, we have used the substitution v = u — t. Since S P; h(i^', 0) = 0, 
this is the desired result, completing the proof. ■ 

As already mentioned, fractional Hougaard motion is our main example of general self- 
similarity outside of the Levy process case. However, the processes S p< h(ij>]s) do not have the 
same marginal distribution as the Hougaard Levy processes of Section 13.21 In particular, for 
p > 1, S P; H(fJ>; s) has support on the whole real line rather than on [0, oo) or (0, oo). This follows 
from the fact that the parameter h is restricted to negative values, which implies that Wh(t,u) 
is negative for u < and positive for < u < t, so that S Pj h{h] s) has support on the whole real 
line. 



5 Lamperti transformations 

We shall now introduce an extension of the Lamperti transformation to the case of general 
self-similarity. Consider a strictly H-SS stochastic process X(t), and recall that the Lamperti 
transformation yields a new strictly stationary stochastic process {Y(t) : t G M} by means of the 
following exponential scaling transformation, 

Y(t) = e- m X{e l ) for t G R, (5.1) 

Lamperti ( 19621 ) and Embrechts and Maeiimal ( 20021 . p. 11). Now, let X be a general H-SS 



see 

family, and define the family of processes Y{p,; t) for \i G O, by 

Y(p; t) = e~ tH X{pie t{H - l) ; e*) for t G R. (5.2) 

Equation (|5.2|) may be obtained from (|2.3|) by taking t = 1 and then substituting c = e . The fact 
that the right-hand side of (|2.3[) does not depend on c suggests that Y(p,; t) might be stationary. 
By the self-similarity of X we obtain, for each \i G O, 

Y(p;t) = e-^XfjJte*^- 1 ^;^) ~ ^(e-^^W^^je-V) = 1) for f G R. (5.3) 

This result shows that the marginals of Y((i;t) are stationary as a function of t for each \i. 

Note that in the special case of a process with drift, X([i;t) = X(t) + fit, say, the process 
(|5.2p takes the form 

t) = e~ tH X{e t ) + fi = Y(t) + \i for t G R, (5.4) 
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which is simply the process (|5.ip with a shift pel. If we assume that the drift term fit is 
added pathwise to X(t), then the shift [i is similarly added pathwise to Y(t). In this case, (|5.2|) 
is equivalent to the conventional Lamperti transformation. In the special case where X is a 
Brownian motion with drift, the corresponding family Y consists of shifted Ornstein-Uhlenbeck 
processes, and (|5.4p represents a class of generalized Ornstein-Uhlenbeck processes. 

Let us now consider the random field {X(/j,;t) :f>0,/i£O}, and assume that (|2.3p holds 
in the sense of equality of the finite-dimensional distributions of the two random fields. In this 
case, the stationarity of the process Y(fi; •) may be shown as follows: 

Y(n s + t) = e-^ H X(^ s+t ^ H -^-e s+t ) 

= e-We-'BXbie'V-VeW-Vie'e*) 

±e- tH X(ne t ( H - 1 1;e t ) 
= Y(n; t) for t G R, 

for all fi G f2. The stationary processes Y(n; •) are generalized Ornstein-Uhlenbeck processes. 

Let us now turn to the inverse Lamperti transformation. Suppose that {Y(fi; t) : \i € 0, t G R} 
is a family of strictly stationary processes. For a given H G R, we define the family X by 

X(fi; t) = t H Y{^it l ' H ; log t) for t > (5.5) 

for all (j, G 0. In the special case where Y(p>; t) = Y{t) + /i, namely a stationary process Y{t) 
plus a shift /j £ 1, then X has the form X{fj,;i) = t H Y(\ogt) + ^t, which is a strictly fZ-SS 
process plus a drift term. 

The discussion of the self-similarity of X is similar to be above discussion for the ordinary 
Lamperti transformation (j5.2|) . Let us assume that {Y(/j,;t) t G M} is a t-stationary 

random field, in the sense that 

Y(j2; t + s) = Y(n; t) for fi G Q and t G R, 

for all s G R, where = here means equality of the finite-dimensional distributions of the two 
random fields. Under this assumption we obtain, for X defined by (|5.5j) . 

c H X(fi;t) = (ct) H Y (fit 1 - 11 ; log t) 

= {ct) H Y^c H ~ l (ct) l - H -\ogt) 

= (ct) H Yd**- 1 (ct) x - H ^logc + logi) 
= Xi/jc 11 ' 1 ;^, 

which shows that the family X is H-SS. 



6 Lamperti-type limit theorems 

We shall now consider some generalizations of Lamperti's limit theorem of ( Lamperti . 19621 ) . 
according to which all self-similar stochastic processes appear as large-sample limits of suitably 
scaled stochastic processes. 
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Lamperti's fundamental limit theorem says that if we are given a stochastic process X(t) 
and a positive measurable function A(c) satisfying 



A(c) — > oo as c — > oo, (6.1) 

and 

A(c)~ 1 X(ct) 4- Y(t) for t > 0, (6.2) 

as c — > oo, where the stochastic process Y(t) has non-degenerate marginal distributions for all 
t > 0, then Y(t) is H-SS for some H > 0, and all H-SS processes appear as limits of this form. 
Moreover, the function A is regularly varying with index H, that is, 

A(c) = c H L{c), (6.3) 



where L is a slowly varying function. See e.g. lEmbrechts and Maeiimal (|2002l . p. 14) for a proof. 

Let us first discuss an infinitely divisible asymptotic version of Lamperti's limit theorem. 
Hence, let us replace the condition (|6.ip by 

A(c) — > oo as c i 0, 

and let us assume that (|6.2p now holds for c | 0. Then a straightforward modification of 
Lamperti's proof shows that Y(t) is again H-SS for some H > 0, and f|6.3[) holds with L slowly 
varying at zero. 

Let us consider the processes with drift corresponding to X(t) and Y(t), namely X(fi;t) = 
X(t) + lit and Y(fj,; t) = Y(t) + fit, say. Then (|6.2p implies that 

^(^"^(/zA^c" 1 ;^) = Aic^Xict) + ^ 

A Y"(t) +fit = Y(fi; t) for f > 0, 

as c — > oo or c J, 0, where ^(/u; i) is a general self-similar family of stochastic processes. 

These results motivate the following question. Suppose we are given families of stochastic 
processes X(/j,;t) and Y(fj,;t), and a function A(c) such that 

A(c) _1 X(/i^(c)c -1 ; ci) 4 y (//; t) for t > 0, as c -> oo or c | 0. (6.4) 

Under which conditions does this imply that the family of stochastic processes Y is self-similar? 
We shall not answer this question in full generality here, but clearly (|6.4p includes as special 
cases the different cases of convergence discussed so far. The following partial result shows a 
case where (|6.4p applies to non-Levy processes. 

Let us replace A(c) by c H in (|6.4|) . which yields the condition 

c~ H X(fic H - l ;ct) A Y(jji) t) for t > 0, as c -> oo or c 1 0, (6.5) 



which corresponds to convergence of the action of the renormalization group (|2.5p to the fixed 
point (|2.3p . We then have the fo llowing Lamperti-type limit the orem, based on an adaptation 
of the original proof, as given by Embrechts and Maeiima ( 20021 . p. 14). 
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Theorem 6.1 Consider a family of stochastic processes Y(p,;t). 

1. If there exists a family of stochastic processes X(/j,;t) such that 116. 5|) holds, then the process 
Y(/j,; t) is H-SS for each \i € ft. 

2. If Y(fi; t) is H-SS then there exists a family X(/i; t) satisfying \6.5\) . 

Proof: We firs show Item 1. in the case c — > oo. The proof in the case c \. is similar. By 
(|6.5p . for any t\,...,% > and for continuity points x%, . . . , x^ of {Y(fi; tj), 1 < j < k} and 
{Y{nb H - 1 ; btj), 1 < j < k}, where b > 0, we obtain 

Hm Pr {c~ H X(fic H ~ 1 ;ct j ) < x j: l < j < k} = Pr {Y(n;tj) < Xj,l < j < k} . (6.6) 

Replacing tj by btj, [i by ^b H ~ l and Xj by xjb H we obtain 

= Pr {Y {fib 11 - 1 ; btj) < Xj b H , 1 < j < k} . 

(6.7) 

Since the limits on the left-hand sides of (|6.6p and (|6.7|) have the same value, due to (|6.5|) . it 

follows by comparing the two right-hand sides that Y([i;t) = b~ H Y(fib H ~ 1 ;bt), which is the 
self-similarity of Y. 

Item 2. is trivially obtained by taking X(fj,;t) = Y(/j>;t). This completes the proof. ■ 



lim Pr 

c— >oo 



~ H X i/i {cb) H ' 1 ; (cb) tj\ < Xjb H , 1 



< 7 < k 



7 Exponential variance functions (H = 1) 

Much like in the case of strict self-similarity, the value H = 1 for the Hurst exponent corresponds 
to a degenerate family of processes, at least when second moments are finite. In fact, if the family 
X is SSSI in the sense of Definition 12.11 the scaling (|2.1|) with H = 1 becomes 

X((m; ct) = cX(n; t) for t > 0. (7.1) 

Under second-moment assumptions, this implies that Cov [X(fi; s), X(/j,; t)] = V(fi)st (compare 
with ()2.2ip ). For t = 1 it follows, in turn, that 

Var [X{fi; s) - sX(fx; 1)] = Var [X(fi; a)] + s 2 Var [X(p; 1)] - 2sCov [X(fi; s),X(pi; 1)] 

= V{^)s 2 + V(fi)s 2 - 2sV(fi)s = 0. (7.2) 

This implies that X(/i; •) is a degenerate, straight-line processes X(n; s) = sX(fj,; 1) a.s. for each 
/i G 9,. 

Let us instead propose an extended definition of self-similarity for the case H = 1. We say 
that a family of stochastic processes X is self-similar with Hurst exponent H = 1 (1-SSSI) if 
there exists a function / : M + — > £1 such that for all c > and /j£fi. 

X(/i + /(c); ct) = c [X{n; t) + tf{c)\ for t > 0. (7.3) 

The interpretation of this definition is that a location change for the rate and a simultaneous 
rescaling of time is equivalent to a location and scale change for the process. In the special case 
/(c) =0 we obtain the degenerate case (|7.ip already discussed. 
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Proposition 7.1 Let X be a 1-SSSI family of stochastic processes satisfying |7.3| ) with a non- 
constant f. Assume that X(fi;t) has finite expectation for all /i £ ft and t > 0. If the function 
1 \— > EpT (//,;£)] is continuous on [0, oo), then up to a translation of fi there exist constants a and 
6^0 such that 

E [X(ji; t)] = ae htl + (it for t > and fiEQ. (7.4) 
When a ^ 0, this implies that 

/(c) = ft" 1 log c /or c>0. (7.5) 

Proof: By the proof of Proposition 12.11 the stationarity of the increments for X implies that 

E[X(n;t)) = a(jj) + b(p)t, 

for suitable functions a and b. In view of (|7.3p . this implies that the functions a, b and / satisfy 

a(/i + /(c)) + 6(/i + f(c))ct = ca(ji) + [b(n) + /(c)] 

for c,t > 0. By comparing the left- and right-hand linear functions of t we obtain the equations 

o(m + /(c)) = co(At) (7.6) 
&(M + /(c)) = %) + /(c). (7.7) 

Since / is not constant, the solution to (|7.6p is a(fj,) = ae bfl , say, where a = a(0) and 6^0 are 
constants. When a ^ 0, this implies /(c) = 6~ 1 logc. From (|7.7p we obtain b{fj) = 6(0) + /U, 
where we may take 6(0) = up to a translation of ji. m 

When the family X satisfies the conditions of Proposition 17. 1| we define the centered family 
Xq by analogy with (|2.8|) . that is, 

X ( f i;t)=X( f i;t)-ae b ». 

It is easy to check that the family Xq is again self-similar, now with a = 0. The mean function 
(|7,4p for Xq then takes the form 

E [X (fi; t)] = pit for t > and /x e O. (7.8) 

Unless otherwise stated, we always assume that the 1-SSSI family X is centered. 

Let us now derive the covariance function of an 1-SSSI family X, assuming that the rate \i 
and variance function V satisfy (|7.8p and (|2.10p . respectively. A simple rearrangement of (|7.3p 
with c = i -1 yields 

X(n; t) & tX(jj. + / {t- 1 )); 1) - tf (t- 1 ) for i > 0, 
which gives the following expression for the variance of the process, 

Var [X(fx;t)\ = + / (t" 1 )) = Vifo (7.9) 
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say. The covariance function now looks similar to the general case, namely 



Cov [Xfa a),X(n;t)] = \ [V^s) + V^t) - V^- \t - s\)} , (7.10) 

for s, t > 0. 

Following |j0r genseni (jl997l . Ch. 4), we consider the family of infinitely divisible exponential 



dispersion models Tw K) (/i, b, t), indexed by the parameter b E M and with domain Q = R for 
the rate parameter ju. For given b, the model Tw oa (fi,b,t) is defined by the variance function 
V(/i) = cr 2 e^ for fi E R and a 2 > 0. The case 6 = corresponds to the normal distribution with 
constant variance function, but for b ^ 0, we may consider Twoo(/j, b, t) to be a Tweedie model 
with power parameter p = oo. For b > (b < 0), the model Twoo(oo,b,t) (T'w 00 (— oo, 6, i)) is 
an extreme 1-stable distribution, and we hence obtain the extended domain 0,^ = (—00,00] for 
b > and = [—00, 00) for b < 0. 



The model Twoo(/i, 6, i) satisfies two separate transformation properties (|j0rgensenl . 119971 . 
Ch. 4), where the first is a translation property, 

c + TwooG^M) = Tw 00 (c + /i,6,te fec ) for c E R, (7.11) 

and the second is a scaling property, 

cTwoo(/i, b, t) = Tw oc (c/i, 6c _1 , tc~ 2 ) for c > 0. (7-12) 

The latter leads us to define Soo(fi,b;t) as the Levy process with marginal distribution 

5oo(/x, b;t) ~ tTwoo(^,6,t) 

= Tw^/ii^i -1 ,* -1 ) for t > 0. (7.13) 

Depending on the sign of b, we find that 6*00(00, b; t) (b > 0), respectively S' 00 (— 00, 6; i) (6 < 0), 
are extreme 1-stable Levy processes. The exponential family of stochastic processes Soo(n,b;t) 
for fi E fib may hence be generated from the respective processes S^ico, 6; t) by exponential 
tilting. As is the case for ordinary Hougaard Levy processes, the families S' 0O (//, 6;i) may be 
characterized by self-similarity, as we shall now see. 

Theorem 7.2 Let X be a non- degenerate exponential family of Levy processes. Then X is self- 
similar with Hurst exponent H = 1 in the sense J7.5| ) if and only if X is a Hougaard family of 
Levy processes Soo(n, b; t) with b ^ 0. 



Proof: Let us first show that the family of process Soo(n, b; t) for fi E is 1-SSSI and satisfies 
the relation 

5oo(^ + 6 _1 logc,6;ct) = c {S^^, 6; t) + i&~ 1 logc} for t > 0, (7.14) 

corresponding to (|7.3p and (|7.5[) (with f{c) = 6 _1 logc). Since we have a Levy process, it is 
enough to consider the marginal distribution of Soo(n,b;t). Using (|7.13p and (|7.1ip . we find 
that the left-hand side of (|7.14p has marginal distribution 

Sco(v + b~ l log c, b- ct) ~ Tw 00 (ct(/i + b- 1 logc), b (cty 1 , (cty 1 ). 
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The right-hand side of J7Ii| follows from ([7TT3D . (17111) and (pU2D . 

c{5oo(/x,fe;t) +t6 _1 logc} ~ cTwoo(^t + tfr 1 log c, &t~\ t _1 c) 

= Tw 0Q (ct (/i + ft^logc) ^(ct) -1 ,^) -1 ) 

This implies the 1-SS property (|7.14p . as desired. 

To show that self-similarity with H = 1 characterizes the Hougaard Levy process, let us 
consider an exponential family of Levy processes X with variance function V, such that (|7.3p 
is satisfied. We know that V is positive and analytic on CI, due to the exponential family 
assumption, so / is not identically zero. Since we are dealing with a Levy process, the variance 
of the process is VarX(,u; t) = tV(fx). Comparing this with (|7.9p we obtain the equation tV(fi) = 
t 2 V(fj, + /(i -1 )), or equivalently, with s = i , 

V( f i + f(s)) = sV( t i). (7.15) 

If /(c) = 0, we obtain the degenerate process (|7.ip . which is ruled out by assumption. Hence, 
by the same type of arguments as used in connection with (|7.6p . the solution to (|7.15p is 

V{ji) = V(0)e bfM for p € R (7.16) 

for some 6^0, and consequently f(s) = 6 _1 logs. Taking a 2 = V(0), the variance function 
(|7.16p characterizes the family Tw 00 (/Lt, &, f), and hence the family of Hougaard Levy processes 
Sooi/J-, b;t), completing the proof. ■ 
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